Homogenization of membrane and pillar photonic crystals 
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Abstract 

We study wave propagation and diffraction in a bidimensional photonic crystal with finite height, 
in case where the wavelength is large with respect to the period of the structure. The device is 
made of materials with anisotropic permittivity and permeability tensors. We derive rigorously 
the homogenized system, using the concept of two-scale convergence. The effective permittivity 
and permeability tensors turn out to be that of a two-dimensional photonic crystal with infinite 
height. 
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Photonic crystals, i.e. dielectric or metallic artificial periodic structures, are generally 
thought of as strongly scattering devices, authorizing the existence of photonic band gaps. 
However, their actual electromagnetic behavior when the wevelength is large with respect 
to the period is also interesting, because it can produce strongly anisotropic behaviors, 
plasmon frequencies, or even left-handed materials The study of the properties of 

these structures in this asymptotic regime comes under the theory of homogenization A 
lot of results are by now very well-known both for 2D and 3D structures. In this paper, we 
consider a photonic crystal made of a collection of parallel finite-size fibers embedded in a 
matrix. This covers the case of structures made out of a layer of bulk materials in which 
holes are made periodically (membrane photonic crystal) but also the case of structures made 
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out of nanopillars (pillar photonic crystal |4|, la, |y, |7|, |8|, |9|), or more generally, structures 
composed of fibers with finite length embedded in a matrix. Our point is to derive the 
effective permittivity and permeability tensors of this structure when the ratio between the 
period of the structure and the wavelength of the incident field is very small. We show, 
using the two-scale convergence method, that the effective, or homogenized, permittivity 
and permeability tensors of these structures are the same as that of infinitely long fibers, 
for which we had already derived rigorous results lOj] . For infinitely long fibers, explicit 
formulas can be derived in some cases 0, [ill 12>Il3|. Let us note that our results hold for 
dispersive and lossy materials. 

The set of fibers is contained in a domain Q = O x [-L, L] of M 3 (cf. fig.l). The space 
coordinates are denoted: x = (xi, £2,2:3) and we also denote xj_ = (£i,x 2 )- The period of 
the lattice is denoted by 77 (see fig. 2). We denote by Y the basic two-dimensional cell of 
the lattice. The obstacle in Y is denoted by P. We consider time harmonic fields, the time 
dependence is chosen to be exp(— iut). For a given monochromatic incident field (E*,H*), 
we denote by (E^, H* 7 ) the total electromagnetic field. Our aim is to pass to the limit 77 — > 
and determine the limit of the couple (E v , H v ). In our methodology, we get at the limit a 
true electromagnetic scattering problem, for a given wavelength A and a bounded obstacle Q 
characterized by some permittivity and permeability tensors. This situation is quite different 
from other homogenization techniques, making use of periodization conditions, in which the 
frequency tends to zero, thus not leading to a diffraction problem but rather to an electro- 
static one 3|. Such an approach can sometimes give useful explicit formulas but generally 
leads to complicated formulations. Moreover, it does not handle the boundary effects which 
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in some cases may lead to some miscomprehensions [15]. The relative permittivity tensor 
e v (x) and relative permeability tensor (x) are described by : 

rji ( x ) = e for x G M. 3 \Q j ^ (x) = /i for x G M 3 \fi 

(x) = e° z 3 ) for xj_ G f2 ' j //" (x) = //> x 3 ) for x ± G (1) 

where y — >■ e° (y) = (e?- (y)) and y — > fjP (y) = (y)) are ^-periodic 3x3 matrix 
functions. The domain Q is periodically filled with contracted cells i]Y x [-L, L] (see fig. 
2). 

The total electromagnetic field (E v , H v ) satisfies 

curl = wuquFHP 

(2) 

curl H 5 ? = -iue^W 

and (E v — E l , H'' — H l ) satisfies Silver-Miiller radiation conditions. 

In order to describe this problem, we will rely on a two-scale expansion of the fields. 
That way, the physical problem is described by two variables: a macroscopic one x and a 
microscopic one y representing the rapid variations of the material at the scale of one basic 
cell, measured by r), that is, at the scale of r]Y. By noticing that there are no rapid variations 
in the vertical direction X3, the microscopic variable is set to be: y = xx/77. Differential 
operators with respect to variable y are denoted with a subscript y. The fields are periodic 
with respect to that microscopic variable (this corresponds to the neighborhood of the Y 
point in the first Brillouin zone). The limit problem obtained by letting 77 tend to 0, will then 
depend on the macroscopic, physical, variable x but also on the microscopic, hidden, variable 
y. The total limit fields will read E° (x, y) and H° (x, y) and the observable, physical, fields 
will be given by the mean value over the hidden variable y: E (x) = \Y | _1 j Y E° (x, y) dy 
and H [x) = \Y\~ j Y H° (x, y) dy, where \Y\ is the measure area of Y . In order to lighten 
the notations, we denote by brackets the averaging over Y, hence H (x) = (H°) and E (x) = 
<E°>. 

The main mathematical tool that we use is a mathematically clean version of the mul- 
tiscale expansion, widely used in various areas of physics. More precisely, for a vector field 
¥ v in (L 2 (fi)) 3 , we say, by definition, that F v two-scale converges towards F° if for every 
sufficiently regular function <p (x, y), F-periodic with respect to y, we have: 

/ E" (x) 4 (x, xj_/e) dx -> [[ E° (x, y) 4 (x, y) dxdy, (3) 
Jn J JnxY 



as T] tends to 0. 

The limit field F° is square integrable with respect to both variables x and y and is 
F-periodic in the y variable (it belongs to L 2 (fl; [L\ (F)) 3 j). A complete analysis of this 



new mathematical tool can be found in 
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We make the physically reasonable assumption that the electromagnetic energy remains 
bounded when rj tends to 0, which is equivalent to assume that (E v , H v ) are both locally 
square integrable. Then it is known [16[ that (E v , H^) two-scale converges towards limit 
fields (E°,H°). This physical assumption could be justified mathematically, however it 
seems quite obvious, from the point of view of physics, that the limit fields exist. The point 
is then to give the system of equations that is satisfied by these fields and to derive the 
effective permittivity and permeability tensors. 

First of all, we have to determine the set of equations that are microscopically satisfied, 
that is, satisfied with respect to the hidden variable y, for that will give the constitutive 
relations of the homogenized medium. Multiplying Maxwell-Faraday equation by a regular 
test function ^x, ^ J , and integrating over Q, we obtain: 

/ E" (x) . curl,, ((f)) + -curl y (0) rfx = iufi I (x) W (x) (x, x±/r]) dx. (4) 
Jn L V J Jn 

Multiplying by 77 and letting r\ tend to 0, we get using Q: 

ft E° (x, y) .curl y (0) dxdy = 0. (5) 
J JnxY 

This is equivalent to: 

/ / curlyE (x, y) .0 (x, y) dxrfy = (6) 
J JnxY 

which is nothing else but the variational form for: curl^E = 0. In a very similar 
way, but using now Maxwell-Ampere equation, we get the equation: curl^H = 0. On 
the other hand, since e v E ri is divergence free, we have, for every test function 0(x, y), 
f n e v (x) ~E V (x) . V x + - Vy0 dx = 0. Multiplying by 77 and letting 77 tend to 0, we get: 

// ^(y)E°(x,y).Vy0rfxdy = O, (7) 
J JnxY 

this can be written as (notice that the div y operator acts only on the transverse components): 

divy (e°E°) = 0. (8) 
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Similarly, since the magnetic field is divergence free, we derive: 



div^ (/i°H°) = 0. 



(9) 



Summing up, we have the following microscopic equations, holding inside the basic cell Y: 



div y (/i°H°) = ] div y (e°E°) = 
curLE = 



curl^H = 



(10) 



The systems in (11 Oj) are respectively of electrostatic and magnetostatic types. This means 
that, with respect to the hidden variable y, the electric field and magnetic field satisfy the 
static Maxwell system. This property is true only at that scale and not at the macroscopic 
scale. However, it is these static equations that determine the effective permittivity and 
permeability. Indeed let us analyze this system starting with the electric field. From the 
curl relation, we get V y E^ = 0, and so E®(x., y) = £"3 (x). Besides, the basic cell having the 
geometry of a torus, we get the existence of a regular periodic function we (y) such that: 



E° = E ± + VyWE- 



(11) 



The function we is the electrostatic potential associated with the microscopic electrostatic 
problem. Inserting (fTT|) in equation (jHJ) and projecting on the two horizontal axis, we obtain: 



div y [e° (a, + V y w E ,i)] =0, i G {1,2} 



(12) 



By linearity, denoting Ej_ = (Ei,E 2 ), we derive that the potential we is given by we 
E\We,i + E 2 we,2, where WE,i are the periodic solutions of (fT2*|) . Thus by (fTTj): 



E° (x, y) = £ (y) E (x) 



(13) 



£(y) 



(14) 



where: 

/ 

+ d yi w E ,x d yi w E ,2 

d y2 w E ,i 1 + d y2 w E ,2 

, 1 

The magnetic field H° can be handled in the same way since it satisfies exactly the same 
kind of equations as H° (see (jlOjO . In particular, we may represent its tranversal component 
in the form: = Hi + Vj_wh, where wh is a periodic magnetic potential (the possibility 
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of this representation is due to the curl-free condition which means that no microscopic 
current is present). Analogously as in (jl3H4j) . we find: 



H° (x, y) = M (y) H (x) 



(15) 



with 



M{y) = 



1 + d yi w HA d yi w Hi2 

8y 2 W H ,l 1 + 9y 2 W H ,2 



\ 



V 











(16) 
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where: 

div y [fj? (e t + V y w H ,i)\ =0,i6 {1, 2} (17) 

The above results show that, at the microscopic scale, the limit fields (E°, H°) are completely 
determined by the physical fields (E, H). Now that the microscopic behavior is precised, we 
are able to determine the macroscopic system satisfied by (E, H). To that aim, let us choose 
regular test function <ft (x) independent of variable y. From Maxwell equations we get 



are 
a 



f n H" (x) .curl (0) dx = -icue J n e n (x) E" (x) (x) rfx 
/ n E 5 ? (x) .curl (0) dx = iufiQ f Q (x) W (x) (x) dx 



(18) 



passing to the limit rj — > 0, we get: 

JJ nxY H° (x, y) .curl (0) dxdy = -iue Q JJ QxY e° (y) E° (x, y) 

JLxy E ° ( x ' y) xurl (0) d * d y = iuj vo Hqxy ^ (y) H ° ( x - y) 

Recalling that (E°) = E and that (H°) = H, we get: 

curl E = iuj^o (fi°H°) 
curl H = -iue (e°E°) 



(x) d'x.dy 
(x) dxdy 



(19) 



(20) 



which, taking into account ()1HJ15|) . brings to the limit system: 

curl E = icufio (fJ°M) H 
curl H = -iue (e°£) E 



(21) 



The homogenized permeability and permittivity tensors are thus respectively (fi°Ai) and 
(e°S). It appears here that the homogenization process is purely local, and that the finiteness 
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of the fibers does not play any role in the homogeneous properties of the medium: the 
effective tensors coincide with that obtained in the polarized cases (3, 10, 17]. This surprising 
property is easily foreseen by our methodology. An approach relying on explicit calculations, 
for instance using Bloch-waves theory or Fourier-Bessel expansions, cannot work here, due 
the lack of an explicit representation of the fields in case of finite size fibers. It should also 
be noted that the case of materials with losses is handled by our result. This result can be 
straightforwardly applied to the study of membrane photonic crystal in the long wavelength 
range where phenomena of birefringence and dichroism are obtained 18]. However, we 
emphasize that the locality pointed out, that is, the fact the effective constitutive relations 
are local ones, is lost if we change the scale of the permittivity coefficients in the obstacles. 
In p articular, the results obtained in the case of infinite conductivities in the polarized case 
|l0j | cannot be transposed to the case of fibers with finite length, due to the emergence of 
surprising non local effects which are studied in (2! . Iz^ . We also remark that the situation 
that we handle here is different from that studied in [20] where the small parameter is the 
depth over wavelength, while the period over wavelength ratio is not small, contrarily to our 
situation. In that case, a dependence on the depth is found. In our homogenization result, 
it is clear that the main numerical problem is the solving of the annex problems ()12I17|) 
for they give the effective matrices £ and M.. In certain simple cases, for instance that of 
circular isotropic non magnetic rods and a permittivity constant in each connected region, it 
is possible to find an explicit expression for the effective permittivity (it is in fact a very old 
problem). However, for more complicated geometries, there is a general numerical procedure 
)ased on fictitious sources, that allows to solve both annex problems at a low numerical cost 

3- 
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Figures captions 

Figure 1: Schematics of the photonic crystal 
Figure 2: Schematics of the basic cells. 

(a) Tridimensional basic cell with cylindrical obstacle. 

(b) Bidimensional basic cell Y with 2D obstacle P. 
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Figure 1 

10 




Figure 2 
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